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Abstract: We use holographic techniques to study SU(N C ) super Yang-Mills theory coupled 
to N { <C N c flavours of fundamental matter at finite temperature and baryon density. We 
focus on four dimensions, for which the dual description consists of N { D7-branes in the 
background of N c black D3-branes, but our results apply in other dimensions as well. A 
non-zero chemical potential /i b or baryon number density n h is introduced via a nonvanishing 
worldvolume gauge field on the D7-branes. Ref. [1] identified a first order phase transition 
at zero density associated with 'melting' of the mesons. This extends to a line of phase 
transitions for small n b , which terminates at a critical point at finite n b . Investigation of 
the D7-branes' thermodynamics reveals that (<9/U b /<9n b )*r < in a small region of the phase 
diagram, indicating an instability. We comment on a possible new phase which may appear 
in this region. 
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1. Introduction 

In strongly coupled, large- N c gauge theories with a gravity dual [2, 3], N f <C N c flavours of 
fundamental matter can be described by N { D-brane probes [4] in the appropriate gravitational 
background. At sufficiently high temperatures, the latter contains a black hole [5]. Working 
with N { <C N c flavours ensures that the matter branes only make a small perturbation to this 
background. Then much of the physics can be studied in the probe approximation where the 
gravitational backreaction of these branes is neglected. 1 

This framework was recently used in [1, 7, 8, 9] to study the thermal properties of N { 
flavours of fundamental quarks (and scalars) in SU(N C ) super Yang-Mills theories in diverse 
dimensions. 2 It was shown that a universal, first order phase transition occurs at some crit- 
ical temperature T fun . At low temperatures, the branes sit outside the black hole in what 
was dubbed a 'Minkowski' embedding (see figure 1), and stable meson bound states exist. 
In this phase the meson spectrum exhibits a mass gap and is discrete. Above some critical 
temperature T fun the branes fall through the horizon in what were dubbed 'black hole' em- 
beddings. In this phase the meson spectrum is gapless and continuous. This large- N c , strong 
coupling phase transition is therefore associated with the melting of the mesons. In theories 

x The backreaction can not be ignored in calculating the effect of the fundamental matter on hydrodynamic 
transport coefficients such as the shear viscosity [6]. 
2 Initial studies include [10, 11]. 
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Minkowski embedding Critical embedding Black hole embedding 




Figure 1: Various possible D7-brane embeddings in the black D3-brane geometry for zero baryon 
number density. The temperature increases from left to right. At finite n b , the Minkowski (and 
critical) embeddings are not allowed - see discussion in the text. 

that undergo a confinement /deconfinement phase transition at some temperature T d < T fun , 
mesonic states thus remain bound in the deconfined phase for the range of temperatures 
T d <T< T fun . 

This physics is in qualitative agreement with that of QCD, in which ss and cc states, for 
example, seem to survive the deconfinement phase transition at T d ~ 175 MeV - see [7] for 
a more detailed discussion. It is thus interesting to ask how this physics is modified at finite 
baryon density. In the presence of N { flavours of equal mass, the gauge theory possesses a 
global U (JV f ) ~ SU(N t ) x J7(l) q symmetry. The U(l) q charge counts the net number of quarks, 
i.e., the number of baryons times N c - see appendix A for details. In the gravity description, 
this global symmetry corresponds to the U(N { ) gauge symmetry on the worldvolume of the 
N { D-brane probes. The conserved currents associated to the U(N t ) symmetry of the gauge 
theory are dual to the gauge fields on the D-branes. Thus, the introduction of a chemical 
potential fi h or a non-zero density n b for the baryon number in the gauge theory corresponds 
to turning on the diagonal U(l) C U(N t ) gauge field on the D-branes. 3 

In this paper we study the gauge theory at constant baryon number density n b . We find 

that, for any finite value of the baryon number density, the Minkowski embeddings, i.e., those 

embeddings where the probe brane closes off above the horizon, are physically inconsistent. 

Hence at finite n b , we focus our study on black hole embeddings. Despite this difference with 

the n h = case, the first order phase transition found there continues to exist here for small 

enough a baryon number density. In this case, however, the transition is between two black 

hole embeddings. For a large enough baryon number density, there is no phase transition as a 

function of the temperature. The phase transition ceases to exist at a critical value n b . These 

results are summarised in fig. 2. This phase diagram also shows a shaded region where the 

black hole embeddings are found to be thermodynamically unstable. While the boundary of 

this region shown in the diagram is qualitative, we have found that the unstable region has a 

limited extent to the left of the line of first order phase transitions. Hence the system must 

find a new stable phase, at least, in this small region - see section 3. 

3 This should not be confused with the chemical potential for R-charge (as considered in, e.g., [12, 13]) 
which is dual to internal angular momentum on the 5* 5 in the gravity description. 
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We focus on four-dimensional J\f = 4 su- 
per Yang-Mills coupled to fundamental mat- 
ter, whose dual description consists of N f D7- 
branes in the background of iV c black D3-branes, 
but our results hold in other dimensions. In- 
vestigations of other holographic systems with 
a chemical potential have appeared previously 
in [14, 15]. An overview of the paper is as 
follows: In section 2, we solve for the embed- 
ding of the D7-branes in the black D3-brane 
geometry. Our discussion includes a brief re- 
view of the black hole background, the equa- 
tions of motion determining the embedding, 
and a careful analysis of the required bound- 
ary conditions. In this section, we also discuss 
the effect of finite n b on the critical solution 
and self-similar scaling found at n b = [1, 7]. 
Finally we present the results of numerically 
solving for the embeddings at various values of 
the baryon number density. Section 3 exam- 
ines the thermal properties of the D7-branes, 

including their stability or lack thereof. Section 4 concludes with a discussion of results. Ap- 
pendix A presents some details of the holographic dictionary relating the worldvolume fields 
describing the D7-brane embeddings to their dual operators in the gauge theory. 



Figure 2: Phase diagram: Baryon number n& 
versus temperature T. The line of first order 
phase transitions ends with a critical point at 
(T f * n ,n b ). The phase which we study is in- 
trinsically unstable in the shaded (red) region. 
This plot shows only a small portion of the full 
phase diagram near the critical point. The ori- 
gin of the axes above corresponds to (n b ,T) = 
(0, 0.986 T fun ). 



2. Holographic framework 



2.1 Black D3-branes 

As first proposed in [2], N = 4 super- Yang- Mills (SYM) with gauge group SU(N C ) is 
holographically dual to type IIB string theory on AdS^ x S 5 with N c units of RR five- 
form flux. The dictionary relating the two sides of the duality equates g s = g 2 M /2i: and 
L 4 /£f = 2g 2 M N c = 2A, where L is the AdS curvature scale - for a review, see [3]. In the limit 
of large N c and large A, the string side of the duality reduces to (weakly coupled) classical 
gravity. At a finite temperature, a black hole appears in the supergravity background [5]. 
Following [1, 7], the black hole metric may be written as 



ds 2 



L 



-Cdt 2 + fdxl 



L 2 

+ -j l d Q 2 + Q 2(in 



where 



Hq) = i 



f(g) = i + 



(2.1) 



(2.2) 



4 This metric is related to the standard presentation with the coordinate transformation g 2 
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The gauge theory temperature is then equivalent to the Hawking temperature of the black 
hole horizon, determined as usual by the surface gravity T = k/2tt. Alternatively, the latter 
temperature may be determined by demanding regularity of the Euclidean section obtained 
through the Wick rotation t — > it E . Then t E must be periodically identified with a period (3, 
where 

This holographic framework allows the thermal behaviour of the strongly coupled gauge theory 
to be further studied with standard semiclassical gravity techniques [16]. In particular, the 
entropy density can be calculated as the geometric Hawking-Bekenstein entropy of the horizon 
[5, 17] 5 

S=-^— = — — = — N 2 T 3 , (2.4) 
AGV X 4G /? 3 2 c ' y ' 

where we have used 167rG = (27r) 7 ^f g 2 . The parametric dependence S oc -/V c 2 reflects the fact 
that the gauge theory is deconfined. Remarkably, this strong coupling result differs from that 
calculated at weak coupling by merely a factor of 3/4 [17]. 

2.2 D7-brane embeddings 

One feature of the JV = 4 SYM theory appearing in the duality above is that all of the fields 
transform in the adjoint representation of the SU(N C ) gauge group. Fields transforming in 
the fundamental representation can be included by introducing an additional set of D-branes 
on the string theory side of the duality. Following [4] , we consider the decoupling limit of the 
intersection of N c D3-branes and N t D7-branes as described by the following array: 

0123456789 
D3: x x x x (2.5) 
D7: xxxxxxxx 

The resulting dual gauge theory is J\f = 4 super- Yang- Mills coupled to N f M = 2 fundamental 
hypermultiplets [4] at temperature T in 3 + 1 dimensions. Assuming N { <C iV c , the decoupling 
limit leads to N { probe D7-branes in the previous background (2.1), with the intersection 
being parametrised by the coordinates {t,x 1 }. Since the D7-branes span the 4567-directions, 
it is useful to introduce spherical coordinates {r, f^} in this space and polar coordinates 
{R, (f)} in the 89-directions. Denoting by 6 the angle between these two spaces we then have: 

g 2 =r 2 + R 2 , r = gsm6, R = gcos9, (2.6) 

and 

dg 2 + g 2 dnl = dg 2 + g 2 (d9 2 + sin 2 6 dQj + cos 2 9 # 2 ) (2.7) 

= dr 2 + r 2 d^l + dR 2 + R 2 d(f) 2 . (2.8) 



5 We divide out by the (formally infinite) three-dimensional volume V x of the Minkowski space in which the 
gauge theory is formulated to yield the (finite) entropy density (2.4). 
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The analysis is simplified by taking x = cos 6 to describe the embedding of the D7-branes. 
Translational symmetry in the 0123-space and rotational symmetry in the 4567-directions 
motivate us to take x = x{o)- The induced metric on the D7-branes is then: 



ds 2 



-Cdt 2 + fdx 2 



+ 



L 2 



l-X 2 + Q 2 (d eX f 

i-x 2 



de 2 + L\i- x 2 )dnl (2.9) 



We also introduce a !7(1) gauge field on the world volume of the D7-branes. As we discuss 
in detail in appendix A, in order to study the gauge theory at finite chemical potential or 
baryon number density, it suffices to turn on the time component of the gauge field, At- 
Again, symmetry considerations lead us to take the ansatz A t = A t (g). The action of the 
D7-branes then becomes: 



-NfTn 




X 2 + B 2 (d e X? - 2(2< 2 ) 2 -^(l - x 2 )F 2 t , (2.10) 



where F gt = d e A t is a radial electric field. 

The equation of motion for A t (Gauss' law) gives 



6 3 f 



(1 - X 2 ?d e A t 



f ^l~X 2 + Q 2 (d e x) 2 ~ 2(2< 2 ) 2 f (1 - X 2 )(d s A t y 



(2.11) 



In the limit that g 
solution approaches 



oo, this equation reduces to d e (g 3 d e A t ) ~ and so the asymptotic 



A t 



a 

fi ~ — + 

g 2 



(2.12) 



The constants \i and a are (proportional to) the chemical potential for and the vacuum 
expectation value of the baryon number density, respectively (see appendix A). The equation 
of motion (2.11) clearly indicates that there is a constant of motion, which we write as 



d = N t T m (2irt\ 



2\2 



Q 3 P 



(l-X 2 ) 2 d e A 



2 / 



(2.13) 



X 2 + Q 2 {d Q x) 2 ~ 2(27r£ 2 ) 2 ^(l - x 2 ){d e A t y 

With this normalization, this constant is precisely the electric displacement, d = SI D7 /8F gt . 
Taking the large-^ limit of eq. (2.13) with the asymptotic form (2.12), we find: 

d = Af f T D7 (2< 2 ) 2 a. (2.14) 

Now one could proceed to derive the equation of motion for the D7-brane profile x{o) 
from the action (2.10) and then use eq. (2.13) to eliminate A t in favor of the constant d. 
Instead, we first construct the Legendre transform of eq. (2.10) with respect to d to eliminate 
At directly at the level of the action. The result is: 

51 



d 8 aF, 



= -N ( T DT J d s a^-ff(l-x 2 )^l-X 2 + Q 2 (d g x) 2 



1 + 



8d 2 



(2.15) 

1/2 



(27rf 2 iV f T D7 ) V/ 3 (l - X 



2\3 
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The gauge field equations resulting from this Legendre transform are simply d e d = 5I D7 /5At 
and d e A t = —SI n7 /Sd. The first of these reproduces the fact that d is a fixed constant and 
we will return to the second one below. 

Before deriving the equation of motion for the D7-brane profile x(f?), it is convenient to 
introduce dimensionless quantities: 

P= i' J= 2^4iV f T D7 - (2 - 16) 
The x equation from eq. (2.15) can then be written as 

(2.17) 



p 5 ff(l ~ X 2 )X L 8d 2 

v±-x 2 +p 2 x 2 v P 6 pa-x 2 ) 

p 3 fh /r~ 



1+ 



3(1 \"! • 2,; 2 \- 2\d 2 1 " ^ + f/ ^ 



Vi-x 2 + p 2 x 2 \ P 6 / 3 (l - X 2 ) 3 L ' p & P{l - X 2 ? + 8d~2j ' 

where the dot denotes derivatives with respect to p, i.e., x = d p x- With p — > oo, this equation 
becomes at leading order: d p (p 5 x) — — 3p 3 x- Hence asymptotically the profile behaves as 

X = - + 4 + -> ( 2 - 18 ) 
P P 4 

where the dimensionless constants m and c are proportional to the quark mass and condensate, 
respectively [1, 7]. The precise relations are given in appendix A. 

Returning to the gauge field, we begin by introducing a convenient dimensionless potential 
and chemical potential: 

A t = ^A t , A = ^V (2-19) 
Then as described above, (2.15) yields the following equation 

d p A t = 2d — J V X <= X . (2.20) 



7(i-x 2 )[p 6 / 3 (i-x 2 ) 3 + 8d2] 

Integrating yields the potential difference between two radii, 



A t (p) - A t (po) = 2d dp . x2 + P 2 * 2 . (2.21) 



' p0 ^f(l-X 2 )[p 6 f 3 (l-X 2 ) 3 + 8<P] 

We will see below that all embeddings of interest extend down to the horizon at p = 1, so 
po = 1 provides a convenient reference point. Further we set At(p = 1) = by the following 
argument: The event horizon of the background (2.1) can be characterized as a Killing horizon, 
which implies that it contains the bifurcation surface where the Killing vector dt vanishes [18]. 
If the potential A as a one-form is to be well defined, then A t must vanish there. Hence we 
can use (2.21) to calculate the chemical potential, i.e., A t (oo), as 

A = 2d[ dp N^-x 2 + p 2 x 2 _ (222) 



V/(l-X 2 )[p 6 / 3 (l-X 2 ) 3 + 8d 2 ] 
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2.3 Near-horizon embeddings 

An important role in [1, 7] was played by the analysis of the probe brane embeddings in the 
near-horizon region of the geometry (2.1). In this section we will see how this analysis is 
affected by the presence of the electric field on the D7-branes. In fact we will generalize the 
analysis to consider probe Dq-branes in a black Dp-brane background, along the lines of [1]. 
These calculations will lead to two main conclusions. The first one is that smooth Minkowski 
embeddings are unphysical for any non-zero baryon density. The second one is that we expect 
the first order phase transition found in [1, 7] to persist for small values of the baryon density, 
but to disappear for sufficiently large densities. 

In order to focus on the near-horizon region, we set 

g = u + —z, 0=f, (2.23) 

and expand the metric (2.1) to lowest order in z, y. This yields Rindler space together with 
some spectator directions which we omit since they will play no role in the following: 

ds 2 = -(2vrT) Vdt 2 + dz 2 + dy 2 + y 2 dti 2 n + ■■■ . (2.24) 

We recall that T = u /itL 2 . In (2.24) we have introduced an integer n equal to the dimension 
of the internal sphere wrapped by the probe Dq-branes. For the D3/D7 system n = 3, but as 
stated above our analysis in this section will apply to more general Dp/Dq systems, for which 
possibly n ^ 3; for example, n = 2 for the D4/D6 system of [11]. The horizon is of course at 
z = 0. The coordinates z and y are the near-horizon analogues of the global coordinates R 
and r in (2.8), respectively. 

In order to describe the embedding of the Dq-branes, we choose the static gauge for 
all their coordinates except the radial coordinate on the brane, which we denote as a. The 
Dq-brane embedding may then be described parametrically as: z = z(a), y = y(cr). We 
modify the analysis of [1] by adding a radial electric field E = £ 2 A t /T, where the dot denotes 
differentiation with respect to a. For simplicity, in this section we will ignore the overall 
normalisation of the Dq-branes action and take I D7 oc f da£, where 

C = -y n y/z 2 {z 2 + y 2 )-E 2 . (2.25) 

This action is homogeneous of degree 2 + n under the rescaling 

z^az, y^ay, E — > a 2 E , (2.26) 

which means that the equations of motion will be invariant under such a transformation. 
Recall that as first described in [19], this scaling symmetry was a key ingredient for self- 
similarity of the brane embeddings in [1, 7]. However, in the present case with E ^ 0, the 
symmetry does not act within the family of embedding solutions with a fixed electric field 
(or rather fixed d - see eq. (2.28) below). Hence we can not expect to find exactly the same 
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self-similar behaviour for branes supporting a fixed chemical potential or baryon density. 
However, we argue below that the embeddings should behave in approximately the same way 
at least where the gauge field is a small perturbation on the Dq-brane. 

As in the previous subsection, it is convenient to work with the electric displacement 

dC y n E 

d = ^= , (2.27) 

8E z 2 {z 2 + y 2 ) — E 2 K 1 

which is constant by virtue of Gauss' law. This is the near-horizon analogue of the quantity 
with the same name introduced in the previous subsection. 6 Note that under the scaling 
(2.26) d transforms as 

d^a n d . (2.28) 
Inverting the relation (2.27) above, one finds 

d 2 + y 2n y ' 



It is also useful to note the relation 



y/z 2 (z 2 + y 2 ) -E 2 = y n z 



z 2 +y 2 ^ l ' 2 



d 2 + y 2n 



(2.30) 



To eliminate E in favour of d and obtain a functional for y(a) and z(a), we perform a Legendre 
transformation by defining 

£ = C- Ed= -z^z 2 + y 2 ^d 2 + y 2n , (2.31) 

in analogy with (2.15). It is easily verified that the equations of motion obtained from C are 
the same as those obtained by first varying £ and then using eq. (2.29) to eliminate E. 

We can conclude from eq. (2.29) that Minkowski embeddings which close off smoothly at 
the y-axis, such as those considered in [1, 7], are unphysical if d ^ 0. These embeddings are 
most appropriately described in the gauge y = a, and they are characterised by the condition 
that the brane reaches y = at some finite z = zq > 0. For the brane geometry to be smooth 
there, we must impose the boundary condition i(0) = 0. Eq. (2.29) then yields E 2 = Zq at 
y = 0. Now even though E remains finite, the tensor field Edy A dt is ill-defined at the origin 
and so one should conclude that these configurations are singular. This singularity is made 
clearer by considering the electric displacement d which also remains constant at the origin. 
However, one should note that as defined in eq. (2.27) d is actually a tensor density and so 



the norm of the associated tensor field is 



d d 8 



= -q^-qI d 2 /y 2n , which clearly diverges at the 
origin. The physical reason why Minkowski embeddings are inconsistent is, of course, that the 
radial electric field lines have nowhere to end if the brane closes off above the horizon. This 
makes it clear that, although we have obtained this result in the near- horizon approximation, 
the same conclusion follows from an analysis in the full geometry (2.1). 



6 Note, however, that they differ in their normalisation. 



- 8 - 



For D-branes, an electric field on the worldvolume can also be associated with fundamen- 
tal strings 'dissolved' into the the D7-brane [20] - see also the discussion around eq. (A. 3). 
Hence the above statement that the electric field lines have nowhere to end can also be viewed 
as the fact that the strings have nowhere to end if the brane closes off. However, rather than 
simply viewing the Minkowski embeddings as unphysical, this point of view lends itself to 
the interpretation that these embeddings by themselves are incomplete. That is, one could 
imagine constructing a physical configuration by attaching a bundle of fundamental strings 
to the brane at y = and letting these stretch down to the horizon. The strings resolve 
the singularity in the electric field since they act as point charges which are the source of 
this field. However, in such a configuration, the strings and the brane must satisfy a force 
balance equation at the point where they are connected. It is clear that if the brane closes 
off smoothly with z(0) = 0, then they can not exert any vertical force in the z direction to 
balance the tension of the strings and so this can not be an equilibrium configuration. One 
might then consider 'cuspy' configurations which close off with a finite i(0) but still at some 
z = zq > 0. In this case, the branes exert a vertical force and so one must examine the con- 
figuration in more detail to determine if the two forces can precisely balance. This analysis 
requires a more careful treatment of the normalisation of the brane action and the fields than 
we have presented here. Hence we defer the detailed calculations to the next subsection where 
we will examine the D7-branes in more detail. However, let us state the conclusion here: no 
Minkowski embeddings can achieve an equilibrium for any (finite) value of i(0). Therefore 
we discard Minkowski embeddings for the rest of our analysis in the following. 

Hence we now turn to consider black hole embeddings which intersect the horizon. Since 
these reach the horizon z = at some y = yo they are conveniently described in the gauge 
z = a. The appropriate boundary condition in this case is then y(0) = 0, and the equation 
that follows from £ is 

(y 2n + d 2 ) [zyy + (1 + y 2 )yy] - y 2n (l + y 2 )nz = . (2.32) 

In view of this equation it is clear that we should expect two qualitatively different behaviours 
for solutions with yfi S> d and j/q <C d. In the first case, it is easy to see that y n » d all along 
the solution, and so we effectively recover the equations of motion for d = studied in [1, 7], 
and therefore oscillatory behaviour around a critical solution for large y: 

y^^/nz + t,, £= n [a sin(alogTz) + b cos(alogTz)] , (2.33) 

{Tz)2 

where a, b are determined by yo. As shown in [1], this oscillatory behaviour eventually leads to 
the property that the quark condensate is multi-valued as a function of the quark mass, and 
hence to a first order phase transition (see figure 5 and the discussion in the next subsection) . 
We thus expect a similar transition if yfi 3> d. 

Incidentally, note that, unlike in the case d = 0, here the 'critical solution' y = ^fnz is 
not an exact solution of eq. (2.32) but only an approximate solution for large y. In particular, 
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there is no exact solution of the form y oc z that just touches the horizon except the y = 
solution. Note also that for black hole embeddings eq. (2.29) gives E ~ z as z — > 0, leading 
to a well defined tensor field at the horizon z = 0. 

We now turn to the case yj <C d, for which the equation of motion (2.32) reduces to 



zy + (l + y 2 )y~0, (2.34) 



whose exact solution is 



(2.35) 



y = yo + z 1 log ^ + y/z 2 -zfj c , (2.36) 



where yo an d z\ are integration constants. Recall that the boundary conditions should be 
y(z = 0) = yo and y(z = 0) = 0. It is impossible to satisfy these conditions with the logarithm 
in eq. (2.36). It is also clearly seen in eq. (2.35) that the general solution is problematic (at 
z = z±) unless z\ = 0. Hence the only physical solution in this regime is precisely the constant 
solution: y = yo- 

Further, we note that the embedding starts very near the horizon with y = yo where 
yo <C d and so we ask how it makes a transition to some more interesting profile of the full 
equation (2.32) far from the horizon. The point is that the term ny 2n z will eventually grow 
large and require y to deviate from a constant. Quantitatively, one finds that the transition 
occurs for z ~ yo {d/y^ ) where the leading solution has the form 

Hence we see the 0(z 2 ) correction to the constant embedding is enormously suppressed in 
this regime y$ <C d. Note that at z ~ yo (d/Uo), the second term is comparable to the first 
and so the Taylor series is breaking down. However, at this point, we still have y n <C d and 
y <C 1. In summary, the solution in this regime is a long spike that emanates from the horizon 
almost vertically, resembling a bundle of strings. 

The analysis above thus leads to the following physical picture. If d is small enough, then 
there is a set of embeddings in the near-horizon region for which y$ 3> d, whose physics is 
similar to that of the d = case. In particular, we expect a first order phase transition to occur 
as a function of the temperature. As d increases, the region where the condition y^^>d holds 
gets pushed outside the regime in which the near-horizon analysis is applicable, suggesting 
that the phase transition as a function of temperature should cease to exist for sufficiently large 
d. This is precisely what the phase diagram in figure 2 confirms. In contrast, the condition 
yo <C d can always be met in the near-horizon region, indicating that solutions for which 
the part of the brane near the horizon behaves as a narrow cylinder of almost constant size, 
resembling a bundle of strings, exist for all values of d. This is also confirmed by our numerical 
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analysis in the full geometry (as illustrated in figure 3), since such type of embeddings can 
always be realised, for any fixed d, by increasing the quark mass (or equivalently by decreasing 
the temperature). In the next subsection we analyse some properties of these embeddings 
more closely. 

2.4 Strings from branes 

The near-horizon analysis above revealed the existence of solutions for which the brane re- 
sembles a long narrow cylinder that emanates from the horizon. One's intuition is that this 
spike represents a bundle of strings stretching between the asymptotic brane and the black 
hole. Examples in which fundamental strings attached to a D-brane are described as an elec- 
trically charged spike solution of the DBI action are well known in flat space [21], in AdS 
space [22] and in other brane backgrounds [23] . Here we would like to formalise this intuition 
by investigating the core region of our D7-brane embeddings in more detail. This analysis 
allows us to investigate the boundary conditions for the Minkowski-like embeddings in detail. 
We begin by rewriting the Legendre-transformed action (2.15) as 

Now recall that \ = cos 9 - see eq. (2.7) - and consider the last factor in the integrand. If 
the embedding is very near the axis, i.e., % ~ 1, then the second contribution in this factor 
can be neglected and eq. (2.38) becomes 

I m - -«^^2 / dtdg^Lj-Jl + g^r 

= " n ^^2 / dtdgy/-gtt {g ee + gee(d e e) 2 ), (2.39) 

where we have used the relation (A. 4) between d and the density of strings ra q . We recognize 
the result above as the Nambu-Goto action for a bundle of fundamental strings stretching in 
the g direction but free to bend away from 9 = on the S 5 . It is interesting to note that 
the term that was dropped provides precisely the measure factor associated with the x l and 
S 3 directions in the limit where the d term vanishes (or is small). In this sense then, the 
D7-brane forgets about its extent in those directions. 

Let us consider the boundary conditions for the configurations which reach the axis 9 = 
at some finite g, i.e., for Minkowski-like embeddings. These embeddings would in general have 
a cusp if d e 9 remains finite at 9 = (a smooth embedding would correspond to d e 9 — > oo). 
As discussed in the previous subsection, to produce a potentially physical configuration, we 
would attach a bundle fundamental strings to the tip of the brane (with precisely the density 
n q ). However, to produce a consistent static configuration, there must be a balance between 
the forces exerted by these external strings and the brane along the ^-direction. The effective 
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tension of the branes can be evaluated in many ways, but here we consider the calculation: 



Now if we wish to calculate the effective tension for a bundle of strings smeared out of the 
x*-directions with density n v the same calculation would apply since eq. (2.39) is precisely the 
fundamental string action. However, these strings would lie vertically along the axis and so we 
would evaluate eq. (2.40) with d Q 9 = 0. Hence for a cusp with any non-zero d e 9, the effective 
tension (2.40) is less than that of the vertical strings. Hence none of these Minkowski-like 
embeddings can achieve an equilibrium with the attached strings for any finite value of d g 9. 7 
We might consider these configurations as the initial data in a dynamical context. Then, given 
the results above, we see that the strings will pull the brane down the axis to the horizon - 
a similar discussion appears in a different context in [24]. In any event, we will not consider 
any of these Minkowski-like embeddings in the remainder of our analysis. 

Now let us consider the black hole embeddings that arise from eq. (2.39). In fact, the 
equations resulting from this action were studied as (a special case of) the string configurations 
describing Wilson loops in the AdS/CFT [25]. In general these solutions are loops which begin 
and end at large g. Hence these are inappropriate in the present context. 8 In this context, 
at finite temperature, there is another class of string configurations, namely strings that fall 
straight into the horizon, which display the screening of the quark-antiquark potential. Using 
this experience, we conclude that the only solutions for eq. (2.39) which reach the horizon 
will be the constant configurations 6 = 6q. Hence, as we saw in the near horizon analysis, the 
black hole embeddings near the 9 = axis are long narrow cylinders of constant (angular) 
cross-section. 

One should ask how far out these constant profiles are valid as approximate solutions 
of the full equations derived from eq. (2.38). The approximation that allowed us to derive 
eq. (2.39) required d 1 / 3 3> psin9, assuming p » 1. Hence the constant solutions 9 = 9q 
should remain approximate solutions out to /^transition ~ d 1 / 3 /9q for small 9q 1. Beyond 
this radius we expect the profile should expand out and approach an asymptotically flat brane. 
However, we can push this transition out to an arbitrarily large radius by taking 9$ —> 0. This 
again suggests that with d 7^ 0, there are D7-brane embeddings which reach the horizon no 
matter how far the (asymptotic) brane is from the black hole. We will verify this result with 
numerical investigations of the full solutions for the action (2.38) in the next subsection. 

Our analysis of the static D7-brane profiles near \ ~ 1 have confirmed the idea that the 
embeddings develop a narrow spike that behaves like a bundle of strings stretching between 
the asymptotic brane and the black hole. It is interesting to extend this idea further by inves- 
tigating the dynamical properties of these spikes. As a step in this direction, let us consider 
our framework with the more general ansatz: x((?, i) an d At(g). 9 After a straightforward 

7 The same conclusion applies for the general Dq/Dp-brane configurations discussed in subsection 2.3. 
8 If we use only a portion of these solutions, i.e., the configuration is cut-off before reaching the loop's 
minimum g, the profile describes the cuspy configurations discussed above. 

9 The symmetries of the problem ensure that this ansatz leads to a consistent solution. 
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calculation the Legendre-transformed action becomes 

1/2 



(2.41) 

As above, we restrict our attention to the embeddings when they are very close to the axis 
~ 0. In this regime, the second contribution in the last factor can be neglected and eq. (2.41) 
becomes 

~ n « v *2^ J dtde vJyT^ 1 + e2{ds0? ~ ^ ( ^ )2 • (2 - 42) 

Once again we recognize this result as the Nambu-Goto action for a bundle of fundamental 
strings stretching in the ^-direction with dynamical fluctuations in the ^-direction. Hence we 
are beginning to see that not just the static properties of the spikes, such as the tension, but 
also their dynamical spectrum of perturbations matches that of a collection of strings; similar 
results have been seen for the dynamics of the Blon spikes on branes in asymptotically flat 
spacetime [26]. In this sense we see that, although no fundamental strings are initially man- 
ifest, the D7-brane spectrum still captures the presence of these strings. This is a satisfying 
result since these strings stretching between the horizon and the asymptotic D7-branes are 
dual to the quarks in the field theory, for which we are turning on chemical potential (i. It 
would be interesting to investigate these issues in more detail. 

2.5 Numerical embeddings 

We now return to the detailed analysis of the D7-brane embeddings in the black D3-brane 
background. In general, it is not feasible to solve analytically eq. (2.17), which determines 
the profile x{p)i so we resorted to numerics. We numerically integrated eq. (2.17), specifying 
boundary conditions on the horizon p min = 1: %(1) = xo f° r various < xo < 1 and 
d p x\p=i = 0. In order to compute the constants m, c corresponding to each choice of boundary 
condition at the horizon, we fitted the solutions to the asymptotic form (2.18). Several 
representative D7-brane profiles are depicted in figure 3. In particular, we see explicitly here 
the formation of long narrow spikes reaching down to the horizon as xo approaches 1 (or R 
approaches 1 on the horizon). 

We can make the appearance of these spikes quantitative here by examining how varying 
boundary the condition xo changes the quark mass m - recall that the latter is proportional 
to the distance which the branes reach along the vertical axis of figure 3. Figure 4 shows 
plots of m versus xo for d = 10~ 4 /4 and 1/4. Note that in both cases, as xo — ► lj the 
quark mass is diverging. Hence with d / 0, there are D7-brane embeddings which reach the 
horizon no matter how large the (asymptotic) separation between the brane and the black 
hole becomes. Since m oc MJT as shown in eq. (A. 2), m — ► oo corresponds to T — > for a 
fixed quark mass M q . Hence the previous result is equivalent to saying that the D7-branes 
intersect the horizon for all values of T when d / 0. Contrast this with the d = case, 
where embeddings of the D7-branes which intersect the horizon {i.e., black hole embeddings) 
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only existed above some minimum temperature [1, 7]. At low temperatures the D7-branes 
were described by embeddings which smoothly closed off above the horizon (i.e., Minkowski 
embeddings). For nonzero chemical potential or nonzero baryon density, there are black hole 
embeddings corresponding to all temperatures in the gauge theory. For small temperatures, 
or large quark mass, most of the brane is very far away from the horizon with only a very 
thin long spike extending down to touch the horizon. Far from the black hole, this embedding 
would look very much like a Minkowski embedding in the low temperature phase of d = 0. It 
differs only by the narrow spike going down to touch the horizon. 
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Figure 4: Quark mass m versus boundary condition \o on the horizon for (a) d = 10 4 /4 and (b) 
d= 1/4. 




Figures 5, 6 and 7 illustrate the depen- 
dence of the quark condensate c on the tem- 
perature T. Several such plots of c versus T 
with varying degrees of resolution are given 
in figure 5 for small values of the baryon den- 
sity: d = 0, 10~ 6 /4 and l(T 4 /4. In the first 
two plots, the differences between the curves 
is virtually indiscernable. In particular then, 
they all begin to show the spiralling behaviour 
that was characteristic of the self-similar scal- 
ing discovered for d = [1, 7]. Of course, 
section 2.3 argued that these spirals should 
persist to a certain level at small d. Note that 

in the highest resolution plot (the last one in fig. 5), one sees that for d = 10~ 4 /4 the small 
scale spirals have been eliminated. In any event, the plots in figure 5 explicitly demonstrate 
that, for small baryon density d, the black hole embeddings are mimicking the behaviour of 
both the black hole and Minkowski branches of the theory at d = 0. Hence certain features of 



Figure 3: Profiles of various D7-brane embed- 
dings in the D3-brane background for d = 10~ 4 /4. 
The black circle represents the horizon. 



- 14 - 



the physics will be continuous between the theories with vanishing and non-vanishing baryon 
number density. In particular, the spiralling or rather the multi-valuedness of c indicates 
there will be a first order phase transition and so the 'melting' transition found in [1] persists 
to small values of the baryon density. 




0.7705 0.771 0.7715 0.772 0.771225 0.77123 0.771235 0.77121 0.771245 

T/H T/ft 



Figure 5: Quark condensate c versus temperature T/M for d = 0, 10~ 6 /4 and 10 _4 /4 on the black, 
red and blue curves, respectively. At low resolution, these curves are all nearly identical and display 
a similar spiralling behaviour. 

As d is increased, the self-similar, spiralling behaviour becomes less and less pronounced 
and eventually c becomes a single-valued function of T/M. To the best numerical accuracy 
that we could achieve, the critical value at which the phase transition disappears in d* = 
0.00315. Figure 6 shows the c in the vicinity of the transition around this critical value. For 
d = 0.0031, the curve shows a slight S-shape and so a small first order phase transition would 
still occur. For the critical value d* = 0.00315, the curve is monotonic but with a singular 
slope near the center. In this case, the phase transition would be reduced to second order. 
Finally for d = 0.0032, the curve is monotonic with a finite slope everywhere and so the phase 
transition has disappeared. 

For completeness, we also show the behaviour of the quark condensate at much larger 
values of the baryon density in figure 7. Figure 7a corresponds to d = 1/4 where some 
interesting structure still persists around T/M ~ 1, which was where c shows a minimum 
in figure 5 at smaller densities. Figure 7b corresponds to d = 10, where c has become a 
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Figure 6: Quark condensate c versus temperature T/M near the critical point. The solid black curve 
corresponds to the critical baryon density d* = 0.00315. The dashed curves above (blue) and below 
(red) correspond to d = 0.0031 and 0.0032, respectively. 



monotonically increasing (towards zero) function of T. 




Figure 7: Quark condensate c versus temperature T/M for (a) d = 1/4 and (b) d = 10. 

We integrated (2.22) numerically to solve for the chemical potential. Plots of jl ver- 
sus temperature all show an apparent divergence as T/M — > 0, as illustrated in figure 8a. 
However, this behaviour is misleading as we now explain. As discussed in the previous sub- 
sections, a common feature of the D7-brane embeddings at small temperatures is the long 
narrow spike close to the 9 = axis. This spike dominates eq. (2.22) for small T/M and so 
the latter formula can be simplified to 

i t-uom 
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where we have restored the dimensions of the chemical potential and the radial coordinate. 
Hence in this limit, the chemical potential is essentially given by the quark mass, as one might 
have expected. Hence the divergence in figure 8a arises simply because fl oc fi/T, as shown in 
eq. (A. 10). This spurious behaviour is eliminated by plotting ///M q = ■\/2pi/m, as shown in 
figure 8b. The latter plot exhibits the small temperature limit /i/M q — > 1 for T approaching 
zero, as is implied by eq. (2.43). Note that if one calculates \i in the vicinity of the phase 
transition, it shows a multi-valuedness similar to that shown for the quark condensate above. 
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(a) (b) 
Figure 8: Chemical potential for d = 10 versus temperature displayed as: (a) jj, and (b) /i/M q . 



3. D7-brane thermodynamics: Free energy, entropy and stability 

We now wish to study the thermal properties of the fundamental hypermulitplets at finite 
baryon number. Our holographic framework translates this question into one of investigating 
the thermal contributions of the D7-branes on the gravity side. As usual, we use the standard 
technique [16] of Wick rotating the time direction. The Euclidean time circle of the black D3- 
brane background then becomes the thermal circle in a finite temperature path integral, and 
the leading contribution to the free energy is determined by evaluating the Euclidean action. 
As we are interested in the contributions of the fundamental matter, we only study the action 
of the D7-branes. Although evaluating the bulk brane action leads to a formally divergent 
result, the AdS/CFT correspondence provides a prescription to remove these divergences: 
One introduces a finite-radius UV cut-off and a set of boundary counterterms to renormalise 
the action [27]. This approach for the branes is completely analogous to the same calculations 
which are performed for the gravity background [28]. This holographic renormalisation of the 
D7-brane action was discussed in more detail in refs. [1, 7], which we follow closely in this 
section. 
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We begin by writing the Euclidean action for the D7-branes in terms of dimensionless 
quantities, introduced in section 2.2, as 10 



4 ulk = J d 8 aC E =M J dpp*ff(l-x 2 )^l-X 2 + P 2 X 2 -j>(l-X 2 )Al (3.1) 

where N is the normalisation constant introduced in [1, 7]: 11 

= 2n 2 N f T m 4 = XN C N ( T 3 

AT 32 1 ' ' 

The normalisation factor illustrates the fact that the leading contributions of the fundamental 
matter are proportional to N c N f , in accord with the large- iV c counting rules of the gauge 
theory. Note then that these contributions are subleading to those of the adjoint fields which 
scale as iV c 2 - see, for example, the entropy density in eq. (2.4). 

As commented above, this bulk action (3.1) contains large-/), UV divergences. Fortu- 
nately, however, these are the same as in the absence of the gauge field, and therefore no new 
counterterms are required beyond those derived in refs. [1, 7], which take the form 

^j^- = ~ (/) max 4 - 2m 2 p max 2 + m A - Amc) , (3.3) 

where /5 max is the UV cut-off. The regularised D7-brane action is then I E = 7 bulk + I bound . It 
can most simply be written as: 

^ = G(m) - l - [(p min 2 - m 2 ) 2 - 4mc] , (3.4) 
where G(m) is the integral: 




G(m) = I dp\ p 3 ff(l - X 2 )y 1 - X 2 + P 2 X 2 ~ ^(1 - X 2 )M -p 3 + m 2 p| . (3.5) 

The limit p max — ► oo may now be taken, since this integral converges. 

As usual, we wish to identify the action with a thermodynamic free energy. However, in 
the present case, there are various possibilities depending on the ensemble under consideration, 
i.e., the Gibbs free energy for the grand canonical ensemble with fixed \x and the Helmholtz 
free energy for the canonical ensemble with fixed n b . Now experience with similar calculations 
for charged black holes, e.g., [12], suggests that the Gibbs free energy is given by the Euclidean 



10 For simplicity, we have left At untouched here rather than introducing a Wick rotated potential At E = 
—iAt- As is well-known, such a Euclidean potential would have to be treated as an imaginary field in the 
present context because the chemical potential and particle density must remain real constants - see, e.g., [12]. 

11 Note that this constant does not include the three-volume V x along the gauge theory directions. Rather 
in this section we will divide out these factors everywhere and so all extensive quantities are actually densities 
per unit volume; for example, (3.1) is the Euclidean action density. 
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action while the Helmholtz free energy is associated with the Legendre transform of I E . Since 
we wish to work with fixed charge, we would want to work with the latter. 

In the following, we confirm the above expectations. Using the equations of motion, the 
variation of the action reduces to a boundary term: 



5I E = 



dC E 9C E r 
d X dA t 



Pmax 



Combining this with the variation of the boundary action I hound (3.3) yields 



6I E = -2Nc5m - -^5jjL 



(3.6) 



(3.7) 



where n q was defined in (A. 4). Recalling that m = M /T we see that the natural thermo- 
dynamic variables of the Euclidean action are the temperature T and the chemical potential 
p. Hence we must identify I E = (3W, where W(T, p) is the thermodynamic potential in the 
grand canonical ensemble, namely the Gibbs free energy. 

Since we wish to work at fixed charge density, i.e., in the canonical ensemble, we perform 
a Legendre transformation by defining 



I v. — Ik 



n q p 
T ' 



which of course is a function of the temperature and the charge density: 

6I E = -2N c5m+ ^-5n a . 



(3.8) 



(3.9) 



We thus identify I E = f3F where F(T, n q ) is the Helmholtz free energy. 
The bulk part of I E is of course the Euclidean analogue of (2.15): 



^bulk 



j dpp 3 ff(l-x 2 Wl-X 2 + P 2 X 2 



1 + 



P 6 P(l-X 2 ) 3 
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(3.10) 



Since the divergences of this bulk action are the same as those of the d = case, the analogous 
expression to eq. (3.4) is now 



^ = G{m) - \ [(m 2 - l) 2 - 4mc] , 



(3.11) 



where G(m) is the integral: 



G(m) 



dp 



1/2 



P 3 ff(i - x 2 Wi-x 2 + p 2 x 2 i + 



P e P(i-x 2 ) 3 , 



p 3 + m 2 p 



(3.12) 

In both of these expressions, we have replaced p m %n = 1 since all of the embeddings which we 
consider terminate at the horizon. 
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We evaluated the free energy numerically for various d and representative results are given 
in figures 9 and 11. The behaviour of the action versus temperature in figure 9 for d = 10~ 4 /4 
is nearly identical to that for d = — we refer the interested reader to compare with the plots 
presented in [1, 7]. The results for d = 10™ 4 /4 are typical for small values of d with the classic 
'swallow tail' shape. Of course, the crossing point of the two branches coming in from small 
and large T marks the temperature of the phase transition. By varying d, one can then map 
out the phase diagram shown above in figure 2. A more detailed diagram is shown here in 
figure 10. We see here that the first order phase transition occurs along a segment starting at 
T fun /M = .7658 at d = and ending at the critical point at T* n /M = .7629 and d* = 0.00315. 




— 1 1 1 1 0.764 0.766 (.768 0.77 0.751 
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Figure 9: Legendre transform of the action, I D7 , versus temperature for d = 10 4 /4. The phase 
transition temperature is denoted by the dotted vertical line in the second plot. 




Figure 10: Phase diagram: Baryon density d versus temperature T/M. 

For completeness, we show some representative plots for large values of d in figure 11, 
where there is no crossing and no phase transition. Note that these plots show an apparent 
divergence as T — > but this is a spurious effect in analogy to the discussion of the plots for 
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the chemical potential. This artifact is actually present in all of the free energy plots but the 
width becomes very narrow at small d. 
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Figure 11: Legendre transform of the action, 7 D7 , versus temperature for (a) d = 10 1 /4 and (b) 
d = 10. There is no phase transition for these values of d. 

We now turn to the entropy density. This can be obtained by differentiating the Helmholtz 
free energy density F(T, d) = TI E with respect to T as 



dF 
df 



-ttL 



dF 

du ' 



(3.13) 



where we used the relation uq = ttL 2 T. Following the calculations described in [7], one must 
carefully consider all of the implicit no dependence in (3.11). The only new contribution 
comes here from the appearance of d in (3.12) since from (2.16), we can see that 



dd 
du 
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Gathering all the contributions, the entropy can be expressed as 

S ~ 

— = -4G(m) + 24d 2 H(m) + (m 2 - l) 2 - 6mc. 
Here we have defined the integral 

fV 1 ~X 2 + P 2 X 2 



H(m) 



dp 



P 3 f 2 (l~X 
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(3.14) 



(3.15) 



(3.16) 



Comparing this expression to eq. (2.22), we see that H = jl/2d. Hence we may write the final 
result as 



-= = -4G(m) + 12 dfl + (m 2 - 1) 



6mc . 



(3.17) 



We evaluated the entropy numerically for various d and some typical results are given in figs. 
12 and 13. The behaviour of the action versus temperature in figure 9 for d = 10~ 4 /4 is 
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nearly identical to that for d = 0. In particular, near the phase transition point, the curve is 
multi-valued because there are several embeddings with the same values of d and T/Al. We 
refer the interested reader to compare with the plots presented in [1, 7]. Figure 13 shows the 
behaviour of the entropy for larger values of d beyond the critical point. 




Figure 12: The entropy S/Af versus temperature T/M for d = 10 4 /4. The position of the phase 
transition is marked by the dotted vertical line in the second figure. 
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Figure 13: The entropy S/Af versus temperature T/M for (A) d= 10 1 /4 and (b) d = 10. 

The thermodynamic identity E = F + T S = T(I E + S) allows us to determine the 
contribution of the D7-brane to the energy density: 
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-3G(m) + \2dfi + 
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(3.18) 



While we did calculate E for many values of d, we do not present any plots here as qualitatively 
their behaviour is similar to that in the plots of the entropy. 
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Finally, we turn to the thermodynamic stability of the system. There are various ways to 
write the requirements for the intrinsic stability of our fixed-charge ensemble. We investigated 
stability here with the conditions: 

§>»• £><>■ <»» 

The first one requires that the system be stable against fluctuations in energy and seems to be 
satisfied everywhere. The second constraint for electrical stability is more interesting, as we 
found that it was not satisfied for all d and T. Our investigations of the region of instability 
remain preliminary, but figure 2 roughly illustrates the extent of the unstable zone as the 
shaded (red) region. In particular, the line of the phase transition seems to be part of the 
boundary of the unstable region between T* un and T fun . This would indicate that the black 
hole embeddings do not correctly describe the true ground state in this small region and in 
particular, just below the phase transition. We hope to return to this matter in the future. 
We comment more on the implications of the instability in the discussion section below. 



4. Discussion 



Ref. [1] identified a universal, first order thermal phase transition in holographic Dp/Dq 
systems. This was characterised by a jump of the Dq-branes between a Minkowski embedding 
and a black hole embedding in the background of the black Dp-branes. In the gauge theory 
this transition is associated to the melting of the mesons. 

Here we have shown that Minkowski embeddings become inconsistent at any finite baryon 
(or equivalently, quark) number density. The physical reason is that a non-zero density which 
is dual to a worldvolume electric field translates into a finite number of strings being dissolved 
into the Dq-branes. Hence the brane is not allowed to close off smoothly as the strings cannot 
simply terminate. We considered the possibility of Minkowski-like embeddings where the 
branes close off above the horizon and external fundamental strings are attached at this point 
and extend down to the horizon. However, examining the forces between the cusp in the brane 
embedding and the external strings, one finds that no equilibrium configuration is possible. 
Rather the strings would pull the tip of the brane down to meet the horizon. We note here 
though that this is not the only possibility for a Minkowski-like embedding. One must simply 
attach a source for the strings and one obvious alternative for such source is the baryon vertex 
[22]. In a Dp-brane background, the baryon vertex consists of a D(8-p)-brane wrapping the 
internal S s ~ p . Hence it may be that there is a family of Minkowski-like embeddings, where a 
gas of baryons absorbs the strings dissolved on the probe branes. It would be interesting to 
investigate this possibility further. 

On the other hand, we did find that with any non-zero baryon density n b , black hole 
embeddings where the Dq-branes intersect the horizon exist for all values of the temperature. 
In contrast, such embeddings do not exist below a certain temperature for n b = and the 
system must be described by a Minkowski embedding beyond this point. In any event, we 
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focused here on studying the behaviour of the black hole embeddings at finite n b in the 
specific example of the D3/D7 system. Our results indicate that the physics is essentially 
continuous around n b = 0. The reason is that black hole embeddings with very small n b mimic 
the behaviour of both n b = Minkowski embeddings and n b = black hole embeddings. 
Moreover, the near-horizon analysis strongly suggested that the universal phase transition 
found in [1] should persist for sufficiently small baryon densities, but that it should cease 
to exist above some critical value n h = n*. This was confirmed by our detailed numerical 
analysis for the D3/D7 system. We emphasize though that the transition at small baryon 
density occurs between two black hole embeddings. 

At zero baryon number density, the spectrum on Minkowski embeddings consists of a 
gapped, discrete set of stable mesons (in the large- N c , strong coupling limit), together with 
stable, massive, free constituent quarks [1, 7]. Instead, mesons on black hole embeddings 
have melted and and the spectrum is continuous and gapless. In fact, little evidence of the 
previously stable states remains in this continuous spectrum [29]. In addition, constituent 
quarks are massless. In the presence of a non-zero baryon density, all embeddings are of 
black hole type and hence no strictly stable mesons exist. Note, however, that the decay 
width is very small if the quark mass is very large, or if the meson is very heavy. Indeed, the 
decay width of a meson is proportional to the support of its wave function on the near-axis 
region where the spike attaches to the branes. This region becomes small as the quark mass 
increases. Alternatively, the peak of the meson wave function occurs further and further away 
from the axis as the meson mass increases - which, for fixed quark mass, can be achieved by, 
for example, increasing the meson radial quantum number. We plan to study these issues in 
more detail elsewhere. 

Similarly, it may seem that the free constituent quarks represent a puzzle in this frame- 
work. Recall that the dual gauge theory is deconfined and so free quarks should play a role, 
in particular since we introduce a chemical potential. The analysis at n b = suggests that at 
least at low temperatures a constituent quark is dual to a string extending from the horizon 
to the brane (at large radius). However, at finite n b , our embeddings are all of the black 
hole type and so if we attach such a string to the brane, it will quickly slip away behind 
the horizon. Hence the puzzle is: How do the D7-branes capture the physics of a gas of 
constituent quarks at low temperatures when there are no stable excitations corresponding 
to macroscopic strings? 

Of course, the resolution of this puzzle is provided by the analysis in subsection 2.4. The 
near-horizon analysis of the Dp/Dq system suggested that, for any value of the baryon den- 
sity, there should exist Dq-brane embeddings which resemble closely Minkowski embeddings 
everywhere except for a long thin spike stretching all the way down to the black Dp-branes 
horizon. This was confirmed for the D3/D7 case by our numerical results, which demonstrate 
that such embeddings correspond to large quark masses (or low temperatures). Further, we 
showed that not only do these spikes match the tension of a bundle of fundamental strings, 
but also their dynamics. Hence these spikes provide a brane realisation of the desired gas 
of constituent quarks. Since the fields describing the D7-branes are dual to meson operators 
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(i.e., operators with n q = 0) in the gauge theory, we may say that, in a very precise sense, 
quarks are being built out of mesons here, in the limit of large quark masses. 

In considering the discussion above, one must remember that part of our phase diagram 
2 corresponds to unstable embeddings. In particular, the line of the phase transition seems 
to be part of the boundary of the unstable region. This would indicate that the black hole 
embeddings do not correctly describe the true ground state of the phase immediately below 
the phase transition. Hence while one should not doubt the existence of a phase transition, 
the precise location of the transition can be called into question. Recall however, that for 
small d^O the behaviour of the black hole embeddings matched everywhere the known 
behaviour of the system with d = very closely, as illustrated in fig. 5. Hence we expect that 
the true line of phase transitions must be very close to that indicated in fig. 10 for small d but 
it may deviate to the right at larger values of d. We also reiterate that we are still refining 
our results on the boundary of the unstable region and that fig. 2 only gives a qualitative 
representation beyond T f * n . It may also be that the region below the phase transition line 
very close to T fun is stable. 

The instability arises in the region where 



It would of course be interesting to identify what the stable ground state is in this region. 
One indication comes from the nature of the instability itself. In the region where (4.1) holds 
the free energy F is a concave function of the baryon density, which means that the system 



One way in which this would be realised in the gravity description would be that it becomes 
thermodynamically favourable for the N { D-brane probes to distribute the U(l) q charge un- 
equally among constituent branes, presumably through some mechanism involving the non- 
Abelian nature of their dynamics. This would imply that the flavour symmetry is sponta- 
neously broken in the infrared. Alternatively, such a separation in different-n b phases may 
be realised by going to a spatially inhomogeneous phase where n b varies from point to point. 
The Minkowski-like embeddings carrying gas of baryons may play a role in this regime. We 
would also note that at this point, it is not clear whether or not other phases or embeddings 
will also play a role beyond the region of instability. In particular, we suspect that a new 
phase may appear at very low temperatures. 

In this paper we have concentrated on the phase structure of gauge theories at constant 
temperature and charge density, namely on their description in the canonical ensemble. It 
will be interesting to consider the phase structure of these theories in the grand canonical 
ensemble, i.e., as a function of the temperature and the chemical potential This should be 
particularly interesting in terms of a potential comparison with the phase structure of QCD. 
However, it is important to keep in mind that much of the interesting physics in QCD at 




(4.1) 




\ such that 



(4.2) 
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finite density - see, e.g., [30] - is associated to the fact that baryon number in QCD is only 
carried by fermionic fields (quarks). This leads to the existence of a Fermi surface at finite 
chemical potential. In gauge theories dual to Dp/Dq systems as those considered here, baryon 
number is also carried by scalar fields, and so the physics at finite chemical potential is likely 
to be very different. In particular, a chemical potential for charged scalars acts effectively 
as a negative mass squared. In the case of free massless scalars this leads to an instability. 
The theories considered here, however, contain interaction, quartic terms in the fundamental 
scalars, and so the chemical potential will presumably lead to condensation of the scalars if 
these are sufficiently light. 
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A. Holographic dictionary 

As described in section 2, the D7-brane embeddings are characterized by two nontrivial func- 
tions, x(p) and A t (p). Further, as usual in AdS/CFT-like dualities, the asymptotic behaviour 
of these fields has a direct translation in terms of operators in the dual gauge theory [3]. In 
particular, considering the asymptotic behaviour in eqs. (2.12) and (2.18), the leading term 
corresponds to the non-normalizable mode and its amplitude indicates the coefficient with 
which the operator is added to the microscopic Lagrangian of the field theory. Similarly, the 
subleading term is the normalizable mode and its amplitude is proportional to the vacuum 
expectation value of the operator. In the present case, since we are discussing worldvolume 
fields, the corresponding operators involve fundamental hypermultiplet fields. 

Let us remind the reader that a hypermultiplet consists of two Weyl fermions ip, if) and two 
complex scalars q,q - the quarks of our theory. These are naturally organized so that ip and 
q transform in the fundamental of the SU(N C ) gauge group, while i/j and q transform in the 
antifundamental. Further, with N { flavours (of equal mass), the hypermultiplets transform 
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under a global U(N f ) ~ SU(N { ) x J7(l) q symmetry. The charges of the fields under the 
diagonal L T (l) q are +1 for ijj and q and -1 for tp and g. Hence the i/(l) q charge naturally 
counts the net number of quarks in a given state. As the colour group is SU(N C ), baryons 
are composed of N c quarks and so we would divide by N c for the number of baryons. 

Now the operators dual to x(p) an d At{p) can be determined by considering the interac- 
tions of the open strings on the D3/D7 array (2.5) before the decoupling limit [32], in analogy 
with the closed strings. Such an exercise leads to the following two operators: 

A t O q = + (q ] V t q - (Vtq^q) + i (q (V t qf - V t qq^ , 

X O m = #V + q(M q + y/2$)tf + q\M q + s/2<f>)q + h.c. . (A.l) 

Recall that the global flavour symmetry discussed above is the U(N { ) gauge symmetry of 
the N f D7-brane worldvolume. Hence O q is simply the quark charge density, i.e., the time 
component of the conserved U(l) q current gauged by on the D7-brane. Note that the 
V t indicate covariant time derivatives in the SU(N C ) gauge theory. The operator O m is the 
variation of the mass term in the microscopic Lagrangian, i.e., O m = —du^- Note that <J?, 
one of the adjoint scalars in the M = 4 supermultiplet, as well as M q , appear in the scalar 
terms here after solving for the auxiliary field constraints within the full coupled theory. As 
a check, one can observe that both of these operators have conformal dimension 12 A = 3, 
which matches the standard prescription for the asymptotic powers appearing in eqs. (2.12) 
and (2.18). 

Now we can make the dictionary between the asymptotic coefficients and the dual gauge 
theory parameters precise by realizing that the hypermultiplet states are the ground states 
of the 3-7 and 7-3 strings. Hence in the decoupling limit, these become precisely strings 
stretching between the D7-brane and the horizon of the D3-brane. For example, the quark 
mass is trivially derived for the brane array (2.5) in asymptotically flat space. As this brane 
configuration is supersymmetric at T = 0, this mass persists in the decoupling limit, where 
it is again the energy of a string stretching between the D3- and D7-branes. This gives a 
relation between M q and the parameter m appearing in eq. (2.18). Further this relation is 
inherited by the theory at finite temperature, since setting T^O does not alter the asymptotic 
properties that determine the gauge theory parameters. 13 Using this result, we can formulate 
a variational argument [11] to relate the second coefficient in the asymptotic expanision of x 
to the field theory condensate (ityip — ^ ^ )) - 1A As the details of this analysis can be found 

12 This dimension applies in the UV where the effects of quark mass are negligible and the theory becomes 
conformal. 

13 Note that here we are referring to the bare or current quark mass. The constituent quark mass is certainly 
modified by thermal effects, as calculated in [7, 33]. 

14 Ref. [11] argued that the scalars would not contribute to the condensate and in writing eq. (A. 2), we have 
certainly ignored such scalar terms. However, strong coupling infrared dynamics might generate a condensate 
for the scalars. See [7] for an interesting example where the scalars dominate this expectation value at high 
temperatures. 
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elsewhere [11, 7], 15 we simply present the results: 

M q = m W> = -2 3/ V 3 £ s 2 iV f T D7 u[j c 

= ^Tm = --VXiV f JV c T 3 c. (A.2) 

8 

Given this result, we note that various figures were plotted in terms of T/M = 1/m and hence 
the relevant mass scale in these plots is M = 2M q /y/\. Up to a factor of 2ir, this corresponds 
to the mass gap in the meson spectrum (at d = 0) [34, 35]. 

Now let us turn to the relation between the D7-brane gauge field and the quark-charge 
operator (A.l). Here the asymptotic value of the potential A t (oo) is proportional to the 
coefficient with which the charge density O q enters the microscopic Lagrangian. This operator 
is normalized so that acting on a particular state it yields exactly the net quark density, and 
therefore the corresponding coefficient is precisely the chemical potential n for the quarks. 
Similarly the relevant expectation value is the quark density n q = (O q ). 

Next we provide a precise definition of the particle density on the string side of the 
duality. First, recall that the electric field on the worldvolume can be thought of as arising 
from fundamental strings 'dissolved' into the the D7-brane [20]. The density of these strings 
can be determined from the local charge density for the two-form i?-field. The standard 
convention is that the fundamental string couples to the NS two-form through the worldsheet 
interaction T { J d 2 a B. Hence a string pointing along the x*-axis sources Bu with the charge 
being just the string tension T { = l/(2ir£ 2 ). Further, the one-form gauge invariance of B 
requires that the D7-brane action only involves the combination B + 2ir£ 2 F. Hence we have 

Slur ^ 1 SI U7 . . 

SB ti ~ 2ir£ 2 6F ti ■ 1 • > 

Combining these observations, we first conclude that since the D7-brane carries an electric 
field in the g direction, the worldvolume effectively contains strings stretching along the radial 
direction with a density precisely determined by the electric displacement d = —5I D r/SFt g . 
The minus sign in the last expression means that for positive d the strings are oriented to be 
inward pointing towards the horizon at g = 1. Since the number of strings corresponds pre- 
cisely to the number of quarks in the field theory, the density of quarks is given by integrating 
the string density on the D7-branes over the internal three-sphere: 

dQ 3 d = 2ir 2 d. (A.4) 

While d is not precisely the coefficient of the normalizable mode in eq. (2.12), the two satisfy 
the simple relation given in eq. (2.14). 

As noted above, the non-normalizable mode At(oo) indicates that the charge density 
operator O q enters the microscopic Lagrangian. As we wish to relate this bulk mode to 



15 Note that the factor of N{ in the formulae for (ipip) was overlooked in [1]. Further [11] only considers the 
case Nf = 1. A full analysis appears in [7]. 
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the chemical potential // for the quarks in the microscopic theory, it is natural to frame the 
discussion in terms of the grand canonical ensemble. There the chemical potential enters the 
partition function as 



cxp 



j d 3 xW([3,fi) =^exp -[3 J d 3 x (H -fiO. 



(A.5) 



where a sum over all states is denoted on the right hand side. Of course, W(/3, /i) and 7i are 
the Gibbs free energy and Hamiltonian densities, respectively. We know that [i oc ^(oo) but 
we would like to determine the exact constant of proportionality. Towards this end, we note 
that, as can be seen from eq. (A.5), 

5W 

^- = -<0 q ) = -n q . (A.6) 

To compare to the string description, we turn to the semiclassical analysis of the Euclidean 
supergravity path integral, as described in section 3. The grand canonical ensemble is repre- 
sented by the usual path integral with fixed At (oo) and the on-shell action gives the leading 
contribution to the Gibbs free energy, i.e., I E = (3 W. Hence to compare to eq. (A.6), we need 
to evaluate the change of the on-shell D7-brane action induced by a change of the boundary 
value A t (oo). Given the worldvolume action (2.10), the desired variation is 

/roc s:r 
dQdn 3 5£ E = 2ir 2 ^ dgj^-d e SA t , (A.7) 

where C is the D7-brane Lagrangian density. In eq. (A.7), we have only integrated over the 
internal three-sphere and the radial direction to produce the free energy density in the gauge 
theory directions. Once again we recognize the first factor as d = 5C/5d g At = —5C E /5d e At 
(in the current notation - note that we need to distinguish between the 'Lagrangian' densities 
appearing in the Minkowski (2.10) and Euclidean (3.1) actions), which is a constant on-shell. 
Hence eq. (A.7) reduces to 

SW = -2tt 2 d (<L4 t (oo) - SA t (l)) = - n q <L4 t (oo) (A.8) 

where we used (A. 4) and the fact that A t always vanishes on the horizon so that we must 
have 5At(l) = 0. Finally comparing eqs. (A.6) and (A.8), we find 

A t (oo)=/i, (A.9) 

and so, as anticipated in the main text, the constant part of the asymptotic gauge potential 
is precisely the chemical potential for the quarks. If we wish to express results in terms of a 
baryon chemical potential, we would convert fi h = N c fi. 

Let us also recall the formulae for the dimensionless quantities defined in eqs. (2.16) and 
(2.19) and which appear in our calculations: 

2tt£ 2 u [2 a , . 

" = ^r = VAf- (A - 10) 

J- d _ _ 25/2 25/2 fAn) 

27r^ n 3jV f T D7 N { N C X 1/2 T 3 JV f AV2 T 3 { ' ' 
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Hence as with the previous definitions, the temperature T provides the scale to make these 
quantities dimensionless but implicitly we have also introduced interesting factors of the 't 
Hooft coupling, as well as of N f and N c . In particular, we see that d is naturally related to 
the expectation value of the baryon number n b in (A. 11). 
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